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ABSTRACT 

We report on the discovery of an ultrasoft X-ray transient source, 3XMM J152130.7-1-074916. It 
was serendipitously detected in an XMM-Newton observation on 2000 August 23, and its location 
is consistent with the center of the galaxy SDSS J152130.72+074916.5 {z = 0.17901 and = 866 
Mpc). The high-quality X-ray spectrum can be fitted with a thermal disk with an apparent inner disk 
temperature of 0.17 keV and a rest-frame 0.24-11.8 keV unabsorbed luminosity of ~ 5 x 10^^ erg s“^, 
subject to a fast-moving warm absorber. Short-term variability was also clearly observed, with the 
spectrum being softer at lower flux. The source was covered but not detected in a Chandra observation 
on 2000 April 3, a Swift observation on 2005 September 10, and a second XMM-Newton observation 
on 2014 January 19, implying a large variability (>260) of the X-ray flux. The optical spectrum of 
the candidate host galaxy, taken ~11 yrs after the XMM-Newton detection, shows no sign of nuclear 
activity. This, combined with its transient and ultrasoft properties, leads us to explain the source 
as tidal disruption of a star by the supermassive black hole in the galactic center. We attribute the 
fast-moving warm absorber detected in the first XMM-Newton observation to the super-Eddington 
outflow associated with the event and the short-term variability to a disk instability that caused fast 
change of the inner disk radius at a constant mass accretion rate. 

Subject headings: accretion, accretion disks — galaxies: individual: 3XMM J152130.7+074916 — 
galaxies: nuclei — X-rays: galaxies. 


1. INTRODUCTION 

When stars wander too close to the supermassive black 
holes (SMBHs), which are believed to be present in the 
center of most massive galaxies, they could be disrupted 
and subsequently accreted, leading to en ergetic flares 
(jLidskii fc Ozernol 119791 iReesI [19^ 11990( 1. Such tidal 
disruption events (TDEs) provide a unique way to find 
and study otherwise dormant SMBHs in galactic nu¬ 
clei. TDEs involving solar-type stars are expected to 
rise very fast (less than a few months) and decay for 
months to years, approximately as after a pe¬ 

riod of peak accretion at a super-Eddington rate. The 
spectrum is expected to be ultrasoft, with character¬ 
istic temperatures of < 0.1 keV, thus mostly in UV 
to soft X-rays. The peak X-ray luminosity Lx can 
reach around 10“^^ erg s“^. Most candidat es discovered 
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X-ray luminosity >10^^ erg s“^) were also found re¬ 
cently and attri buted to the presence of a relativistic 
jet in the events (fBloom et al.ll201lUBurrows et al.[2MTI : 


Cenko et al.l I20l3l . We refer to iKomossal ( 2012 ) and 

KomossJ l 20151) for recent reviews of TDEs in X-rays. 


There are other X-ray transients reaching X-ray lumi¬ 
nosities comparable to those of TDEs around SMBHs. 
However, their duration and lightcurve evolution is typi¬ 
cally very different. The X-ray afterglow of 7-ray bursts 
(GRBs) can hav e peak Lx 10 ^^ erg and last 
for months (e. s.. iLevan et al.l 1201411 or, in a few cases, 
for years fe.g.. iGrupe et al.ll2010ll . Their X-ray spectra 
are g enerally hard with photon index <2 ()GruDe et al.l 
l2013fl . Recently, several ultralong GRBs were discovered, 
with the late-time X-ra y spectra found to be very soft 
(phot on index >3, e.g., (Piro et aI1l20l4 iMargutti et al.l 
I2OI5II . While sup ernovae (SNe) seldom have Lx a bove 
10'^^ erg s“^ le.g.. Ilmmleill200'7l : ILevan et al.ll2013|l . the 
superluminous supernava SGP 06F6 was detected in X- 
rays ^^1 50 days after its initial d iscovery, with Lx ~ 10^^ 
erg s“^. iSoderberg et al.l (|2008ll also serendipitously dis¬ 
covered an extremely luminous X-ray outburst from SN 
2008D, with peak Lx 6 x 10"^^ erg s“^. The main 
event, however, last only several hundred seconds and 
could arise from the prompt shock breakout from the SN 
progenitor’s surface. The outburst showed a significant 
hard-to-soft spectral evolution, with photon index ~1.7, 
typical of SNe, at the peak, to ~3.2 about 400 s later. 

In our continuing effort to classify X-ray sources 
seren dipitously detected by XMM-Newton and Chandra 
(e.s.. iLin et I20TI . we discovered a supersoft 


X-ray transient 3XMM J152130.7+074916 (J1521+0749 
hereafter) in the XMM-Newton Serendipitous Source 
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Catalog (the 3XMM-DR5 version, iRosen et al.l 120151) . 
The source seems to show an X-ray outburst in 2000. 
In this paper, we report the properties of this source 
and argue that it is probably a soft TDE at a red- 
shift of z = 0.17901 (the source luminosity distance 
dL = 866 Mpc, assuming a flat universe with Ho=70 
km s“^ Mpc“^ and r2M=0.3). In Section [21 we describe 
the data analysis of XMM-Newton, Chandra, Swift and 
ROSAT observations. In Section |3l we first identify the 
host galaxy of our source, followed by the presentation 
of its detailed X-ray spectral and timing properties. We 
discuss the nature of our source and give conclusions of 
our study in Sectional 

2. DATA ANALYSIS 

Table [1] lists all the pointed X-ray observations that 
covered J1521-I-0749 and were analyzed by us. There are 
two XMM-Newton observations, one Chandra observa¬ 
tion, one Swift observation, and two ROSAT observa¬ 
tions. They are denoted as XI, X2, Cl, SI, R1 and R2, 
respectively, hereafter (refer to the table). All these ob¬ 
servations have the cluster of galaxies MKW 3s as the 
target, with J1521-I-0749 observed at off-axis angles of 
^6'-9'. J1521-I-0749 was serendipitously detected in XI 
but not in the other observations. We extracted the spec¬ 
trum and light curve from XI for detailed study. We also 
extracted the spectra and the corresponding response 
files for the other observations in order to constrain the 
luminosity level of our source in these observations. Cir¬ 
cular regions (see Table [T] for the radii used) were used 
to extract the source spectra. Our source is relatively far 
away from MKW 3s, but there is still some small contam¬ 
ination of the X-ray emission from this cluster of galaxies 
on our source. In order to accurately represent the con¬ 
tamination level, we extracted the background spectra 
from 4-8 circular regions that are near our source and 
had the same size and the same distance to MKW 3s as 
the source region. Below we give more details about the 
procedures that we adopted in analysis of each observa¬ 
tion. 

The source was in the field of view (FOV) of all 
the three European Photon Imaging Cameras (i.e., pn, 
MOSl, and MOS2. iJansen et ai.1 1200 It iStriider et alJ 
120011 [Turner et al.l 1200 ID in the imaging mode in both 
XMM-Newton observations. We used SAS 13.5.0 and 
the calibration files of 2014 August for reprocessing the 
X-ray event files and follow-up analysis. The data in 
strong background flare intervals, seen in all cameras in 
XI and in the pn camera in X2, were excluded follow¬ 
ing the SAS thread for the filtering against high back¬ 
grounds, i.e., excluding all times when the background 
level exceeded the low and steady leve0- The final ex¬ 
posures used are given in Table |T] In the end we used 
73%, 81%, and 81% of the total exposures of pn, MOSl, 
and MOS2 for XI, respectively. The corresponding frac¬ 
tions are 97%, 100%, and 100% for X2, respectively. The 
event selection criteria that we adopted fo llowed the de¬ 
fault values in the pipeline (see Table 5 in IWatson et al.l 
I2009D . J1521-1-0749 was also covered in the Optical Moni¬ 
tor fOM. [Mason et aD200lD in both observations. In XI, 
the UVWl (2910 A), U (3440 A), B (4500 A), and V 

® http://xmm.esac.esa.int/sas/current/documentation/threads 
/EPIC_filterbackground.shtml 


(5430 A) filters were used, while in X2, two more filters, 
UVWS (2120 A) and UVMS (2310 A^ were used. We 
used the SAS task omichain to obtain the photometry 
of our source (Table [2|). 

The Cl observation used the imaging array of the 
AXA F CCD Imaging Spectrometer (ACIS; iBaiitz et al.l 
I1998D . Our source is in the front-illuminated chip 12. 
We reprocessed the data to apply the latest calibration 
(CALDB 4.5.9) using the script chandra_repro in the 
Chandra Interactive Analysis of Observations (CIAO, 
version 4.6) package. We used a source extraction re¬ 
gion enclosing 70% of the point spread function. 

In the SI observ ation, the X-ray telescope (XRT; 
[Burrows et al.l[2005D was operated in Photon Counting 
mode for 9.6 ks, and the UV-Optical Telescope (UVOT; 
[Roming et al.[[2(i?)^ used the UVWl filter (2910 A) for 
8.9 ks. We analyzed the data with FTOOLS 6.16 and 
the latest calibration files (XRT 20140730 and UVOT 
20130118). The X-ray data were reprocessed with the 
task xrtpipeline (version 0.13.1) to update the calibra¬ 
tion. Other than estimating the luminosity level of our 
source in the XRT, we also calculated the UV emission 
level with the task uvotsource, using radii of 5" and 20" 
for the circular source and background regions, respec¬ 
tively. 

The ROSAT observations were also analyzed with 
FTOOLS 6.16. We note that there are three more 
ROSAT pointed observations that covered J1521-I-0749. 
We did not include them in this study because they are 
not deep enough to provide good constraints on the flux 
limits of our source, due to large off-axis angles and/or 
short exposures. J1521-I-0749 was not detected either in 
these observations. 

3. RESULTS 

3.1. The X-ray Source Position and the Host Calaxy 

The position of J1521-I-0749 in the 3XMM-DR5 
catalog is only 0706 away from the Sloa n Digi¬ 
tal Sky Survey (SDSS, [Abazaiian et al.[ [2009D galaxy 
SDSS J152130.72-1-074916.5 (Figure [TJ. However, such 
a small offset could be an artifact to some extent, con¬ 
sidering that this source was among those used for 
the absolute astrometry correction and that it is one 
of the brightest sources in XI. Therefore, we re¬ 
calculated the source position. We first obtained the 
mean positions (weighted by the uncertainties) of the 
point sources in XI, X2 and Cl, after correcting the 
astrometry of X2 and Cl relative to XI using the 
SAS task catcorr. For Cl, we in cluded a system¬ 
atic u ncertainty of 0716 obtained in [Rots fc Budavaril 
(|2011D . The task catcorr also calculated the system¬ 
atic positional uncer tainties from the astrometric rec¬ 
tification procedure ([Rosen et al.l[201^ . Such system¬ 
atic uncertainties were taken into account and prop¬ 
agated. The mean positions, excluding J1521-(-0749 
and those in bright X-ray emission from MKW 3s, 
were then cross-correlated with optical sources to ob¬ 
tain the absolute astrometry correction. For the optical 
sources, we stacked the Canada-Fr ance-Hawaii Telescope 
(CFHT) MegaPrime/MegaCam ([Boulade et al.[ [2003D 
images (aligned to the SDSS astrometry) in the r' band 
taken betwee n 2008 April and 2010 April using MegaPipe 
(|Gwvn[ [2008D . resulting in a total exposure of 3840 s 
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TABLE 1 

The X-ray Observation Log 


Obs. ID Date 

Detector 

OAA 

T 

Tgrc 

Count rate 

Fabs 

Funabs 




(ks) 


(10“^ counts s“^) 

(10'‘'* erg 3“'') 

(1) (2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

XMM-Newton: 

0109930101(X1) 2000-08-23 pn/MOSl/MOS2 

9.2' 

31/41/41 25" 

/2b"/2b" 

80.6 ± 2.7/16.1 ± 1.1/19.5 ± 1.2 3.457o-2“ 

5 28'*'^'^® 

0723801501(X2) 2014-01-19 pn/MOSl/MOS2 

8.6' 

96/117/117 10" 

/10"/10" 

< 0.17/< 0.18/< 0.12 

< 0.01 

< 0.02 

Chandra: 

900(01) 2000-04-03 

ACIS-I2 

7.6' 

57 

5U" 

< 0.06 

< 0.03 

< 0.04 

Swift: 

00035189001(31) 2005-09-10 

XRT 

6.2' 

9.6 

20" 

< 0.36 

< 0.27 

< 0.42 

ROSAT: 

rp800128n00(Rl) 1992-08-15 

PSPCB 

6.2' 

9.6 

30" 

< 0.99 

< 0.12 

< 0.19 

rh800425n00(R2) 1994-08-23 

HRI 

6.2' 

13.4 

10^^ 

< 0.41 

< 0.22 

< 0.33 

Note. — Columns; (1) the observation ID 

with 

our designation 

given in 

parentheses, (2) the observation 

start date, (3) the 


detector, (4) the off-axis angle, (5) the exposures of data used in final analysis, (6) the radius of the source extraction region, (7) 
the net count rate (0.2—2 keV for XMM-Newton observations, 0.3—2 keV for Swift and Chandra observations, and 0.1—2.4 keV for 
ROSAT observations, all in the observer frame), (8) rest-frame 0.24—11.8 keV luminosity, corrected for Galactic absorption but not 
intrinsic absorption, (9) rest-frame 0.24—11.8 keV luminosity, corrected for all neutral and ionized absorption. The spectral shape 
obtained for XI using the MCD model subject to a fast-moving warm absorber was assumed in calculation of the luminosities for 
other observations. All uncertainties and upper bounds are at the 90% confidence level, calculated using the Bayesian approach with 
the CIAO task aprates for all observations except XI. 


TABLE 2 

The XMM-Newton/OM photometry of J1521+0749 


Obs 

UVW2 

UVM2 

UVWl 

u 

B 

V 

XI 

X2 

> 18.5 

> 19.4 

> 20.0 
> 20.4 

> 20.4 

> 20.6 

21.5 ± 0.5 
21.2 ± 0.4 

20.3 ± 0.5 
20.2 ± 0.3 


Note. — The AB magnitudes with Icr uncertainties or the 
3(7 detection limits are given. 
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Fig. 1. — The CFHT/MegaPrime r'-band image around the 
field of J1521-I-0749. The green circle of radius 0758 (i.e., 1.8 
kpc) represents the 99% positional uncertainty of J1521+0749, 
indicating that J1521-h0749 is consistent with the center of the 
galaxy SDSS J152130.72+074916.5. The two closest galaxies to 
J1521+0749 in MKW 3s are SDSS J152131.80+074851.9 and SDSS 
J152134.15+074815.3 (’gall’ and ’gal2’ in the plot, respectively). 
The inset plot in the upper left corner zooms in on the 5 ” x 5 ” re¬ 
gion around SDSS J152130.72+074916.5. The center of the galaxy 
is marked with a cross. 


and a median seeing FWHM of 0"77. The new absolute 
astrometry correction gave the position of J1521+0749 
to be R.A.=15:21:30.75 and Decl.=+7:49:16.7, with the 
99% positional uncertainty of 0758 (i.e., 1.8 kpc at the 
redshift z = 0.17901 of the candidate host galaxy of 
J1521+0749; see below). The relatively small positional 
uncertainty is mostly thanks to the very high quality 
of the XI detection. This new position is 0733 away 
from the center of SDSS J152130.72+074916.5, consis¬ 
tent within the uncertainty (Figure[T]). The number den¬ 
sity of the optical sources that are as bright as or brighter 
than SDSS J152130.72+074916.5 in the r' band within 
10' is 0.00038 per square arcsec. Thus the chance prob¬ 
ability for our X-ray source to be within 0733 from the 
center of SDSS J152130.72+074916.5 is only 0.0001. 

The SDSS took a spectrum of the galaxy on 2011 May 
28, which is shown in Figure [2] (from the SDSS archive). 
Although the spectrum is a little noisy, the typical ab¬ 
sorption features indicating a non-active galaxy can be 
seen. The redshift of the galaxy is z = 0.17901 + 0.00007 
{Dl = 866 Mpc). No clear emission lines, except 
some spikes due to contamination from sky lines, can 
be seen. The 3a upper limit of the flux of [O 111] A5007 
is 1.8 X 10“^^ erg s“^ cm“^, corresponding to a luminos¬ 
ity of 1.6 X 10^® erg s“^ (Galactic extinction corrected). 
Using the bo lometric correct i on fac tors from the [O 111] 
A5007 flux in lLamastra et al.l (I2009D . we obtained the 3a 
upper limit of the bolometric luminosity of the persistent 
nuclear activity to be 1.4 x 10^^ erg s“^. 

In order to examine the host galaxy proper¬ 
ties in greater detail, we fitted the SDSS spec¬ 
trum to multi-component models comp rised of single¬ 
populat ion synthetic spectra, as in iMaksvm et al.l 
(12014^. Using Penalized P ixel Fittin g (PPXF) soft¬ 
ware (iCaPDellari fc Emsellemi 120041 1 and IVazdekis et all 
(|2010f) synthetic spectra spanning a grid of 48 ages up 
to 14 Gyr, with [M/H]={+0.22, 0.00, -0.40, -0.71, 
—1.31, —1.71}, we found that J1521+0749 is consis¬ 
tent with an old (> 10 Gyr), passive stellar popula¬ 
tion. Full-spectrum kinematics implied a stellar disper¬ 
sion O’* 66 km s~^, which is comparable to the limit 
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TABLE 3 

Fitting results of the XI spectrum from J1521+0749 


Models 

MCD-^Ga'^ 

BB+Ca'^ 

Nn.i (10^“ cm-^) 

0.3+1'^ 

0.0+“+ 

^Tmcd/^’Tbb (keV) 

B I or-l-U.UUO 
U.lZO_Q QQg 

^•'^^^ — 0.003 

Nmcd/Nbb 


6291)““ 

Sga (keV) 

0-75l°:“? 

''^-0.03 

o-ga (keV) 

g Q2 + 0.06 

0.00+““® 

Afga (10-'") 

l-Oio.5 

0.8t“;® 

Labs (10^^ erg 

3.48lS:29 

3.20t“;)“ 

Lunabs (10*^ erg s“^)° 

Q cc + 0.58 

0.00_o.24 


Lbol (10^^ erg 

• •^0_Q g3 

4 24+0-26 

C/u' 

536.4/528 

545.9/528 

Models 

edge*MCD 

edge*BB 

Na.i (10^“ cm“^) 

0.0+“® 

0.0+“+ 

LLmcd/LLeb (keV) 

o-i44tS:SS? 

0.10612:^3 

Nmcd/Nbb 

105lf„ 

458^®^ 

L^edge (keV) 

U.50_q Q2 

0"®4_0.03 

Tedge 

u.yo_o.32 


Labs (10^^ erg 

o or + 015 
0.00_o.i5 

3.07j;“;(“ 

Lunabs (10*^ erg s“^)° 

Q /I d” • 21 

3-46Io.13 

Q ■] Q+0.13 

'5-0.12 

Lbol (10^^ erg 

3Q 

^•k^^_0.20 

C/u' 

529.7/529 

552.8/529 

Models 

zxipcPMCD 

zxipcf*BB 

Nh,! (10^*^ cm“^) 

0.0+ + ® 

0.0+“® 

^Tmcd/^Tbb (keV) 

u.ius;_Q 

0.130iS:TO7 

Nmcd/Nbb 

"70+36 
' -^-20 

577t.% 

L7H,axipcf (10^^ cm”^) 

4 7+2"'^ 

' -1.2 


^zxipcf 


1 29+°°^ 

■^zxipcf 


—0 12+°-°^ 

Labs (10^^ erg 

o 4C + O.I6 
'J-^0_o.29 

8 2n+®'^‘^ 

'5-^U-o.i4 

Lunabs (10^^ erg 

5.28l?;““ 

9 94+^-®’^ 

Lbol (10^^ erg 

8.6j:“;“ 

ii.6j:+,“ 

C/u' 

523.8/528 

532.1/528 

Models 

zxipcP'optxaenfq.-v, ^ 

zxipcPoptxagnfKe-rr^ 

Nn,i (10^“ cm-^) 

0.0++® 

0.0+ + ® 

Mbh (-^g) 

1.9+“-J X 10® 

i.4j:o “ X10® 

T'Boi/ f^Edd 

4 2+^-^ 

Q c-1 +0.28 

^•Oi_o ;12 

AfH,s=.ipcf (10^^ cm-^) 

4 7+2-'^ 

' -2.3 

4 7+2-’^ 

' -2.3 

fog ^zxipcf 


^-'-^‘^- 0.07 

•21zxipcf 


—0 12+°-°^ 

C/u” 

523.7/528 

523.7/528 


Note. — All fits used the C statistic. ®'“Ga” is a Gaussian 
emission line; ^rest-frame 0.24—11.8 keV luminosity, corrected for 
Galactic absorption but not intrinsic absorption; ‘^rest-frame 0.24- 
11.8 keV luminosity, corrected for all neutral and ionized absorp¬ 
tion; ‘^the bolometric luminosity based on the total fiux of the ther¬ 
mal component (MCD or BB); ®the C statistic and the degrees of 
freedom; ^for thermal emission from an accretion disk around a non¬ 
rotating Schwarzschild BH; ®for thermal emission from an accretion 
disk around a maximally rotating Kerr BH. 
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Fig. 2.— Top panel: PPXF best-fit modeling of J1521+0749. 
The SDSS rest-frame data (black, solid) are overplot with the best- 
fit PPXF model (red). Points are smoothed with a boxcar function 
over 5 pixels for clarity. Smoothed residuals are also plotted (blue) 
with mean noise (black, dashed). A strong O I A5577 A (observed) 
sky line has been masked out. Important stellar absorption an<l 
AGN diagnostic emission lines are overplot for reference. Note 
that several strong false ‘features” are associated with atmospheric 
OH complexes at long A, particularly near Ha A6563 A. Middle 
panel: As with the top, with the region about [O Hi] expanded 
for clarity and emphasis. Bottom panel: Relative mass fractions of 
different stellar populations with respect to metallicity and age, as 
per the best-fit model the PPXF template synthesis code. Dom¬ 
inant populations are indicated by shading, with darker shading 
indicating a larger fraction of the total stellar mass in the best-fit 
model. 


from SDSS instrumental dispersion a* ^ 69 km s~^. We 
likewise modeled the spectrum wit h STARLIGHT and two 
ar rays of iBruzual fc Charli^ () 20 n.S[ ) synthetic spectra, as 
in IMaksvm et alJ (l2014bD . STARLIGHT inferred a stel¬ 
lar mass M* ~ 8.4 x 10® Mq. Care should be taken 
in interpreting ct*, since in general S/N < 10 for A < 
5400 A(rest-frame). Using [Graham fc Scottl (j2015ll . we 
inferred the central BH mass Mbh 2 x 10^ Mq, which 
has a 1(7 uncertainty of 0.83 dex. The BH mass can also 
be estimated usin g the BH mass versus bulge X-band lu- 
minosity relation ( Grahan]l2007t iMarconi fc Huntll2003[l 
and is < 5 X lO"^ Mq, based on the X-band magnitude 
of rriK = 16.99 = —22.4) of the galaxy from the 

the 3.8m United Kingdom Infra-red Telesc ope f UKIRT) 
Infrared Deep Sky Survey (|Lawrence et al.l[2007ll . These 
BH mass estimates are consistent with each other and are 
a little larger than the estimate from the X-ray spectral 


fit in Section [3)3] (^10®-10® Mq). 

3.2. The Long-term Variability 

Figure |3| plots the long-term evolution of the unab¬ 
sorbed luminosity Lx of J1521-I-0749 in its rest-frame 
0.24-11.8 keV energy band (or the observer-frame 0.2-10 
keV energy band). The luminosities were obtained based 
on the fit to XI with a multicolor disk (MCD) subject to 
a fast-moving warm absorber fSection l3.3l) and adopting 
a distance of = 866 Mpc ISection [3.1|) . The source 
was not detected in observation Cl in 2000 April, with 
Lx < 4.1 X 10^^ (the 90% confidence upper bound; the 
same below; Table [T]). The source was detected in XI 
4.7 months later, with Lx ~ 5.3 x 10^® erg s“^. In SI 
in 2005 September, the source was not detected, with 
Lx < 4.2 X 10^® erg s“^. It was not detected either in 
X2 in 2014 January, with Lx < 2.1 x 10^^ erg s“^. There- 
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MJD-51637 (day) 

Fig. 3. — The long-term unabsorbed luminosity curve in the rest- 
frame 0.24-11.8 energy band. The uncertainties and upper bounds 
are at the 90% confidence level (see Table fTl. The observations are 
noted in the plot. The solid line represents a (t — decline, 

assuming that XI is on it and that the disruption time to is at the 
time of Cl (i.e., 4.7 months before XI). 



0.3 1 5 

Energy (keV) 


Fig. 4.— (top panel) The unfolded spectrum of XI (for clarity, 
only the pn spectrum is shown). The spectrum was fitted with 
the MCD model subject to a fast-moving warm absorber (solid 
line), (middel panel and bottom panel) The fit residuals with the 
MCD model without and with the warm absorber, respectively 
(pn in black, MOSl in red, and MOS2 in green). For clarity, the 
spectrum was rebinned to be above la per bin in the plot. 

fore the source has a long-term variability factor of >260. 
Figure |3] does not include the ROSAT observations. We 


estimated Lx in R1 and R2 to be < 1.9 x 10^^ erg s“^ 
and < 3.3 x 10“^^ erg s“^, thus a factor of >28 and >16 
lower than that in XI, respectively. 

In the UV and optical, the source was detected in the 
B and V filters, but not in other filters, in XI and X2 
(Table [2]). The optical fluxes in these two XMM-Newton 
observations are consistent with each other and the SDSS 
measurements {u' = 22.4 ± 0.4 mag, g' = 21.09 ± 0.05 
mag, r' = 19.96 ± 0.02 mag, i' — 19.51 ± 0.02 mag, and 
z' = 19.17 ± 0.06 mag). In the UVWl filter in SI, there 
is no clear UV emission from J1521-I-0749 either, except 
some contamination from the read-out streak of a bright 
source. Although we obtained the UVWl AB magnitude 
of 23.3±0.3 mag for J1521-I-0749 using uvotsource, this 
value should be viewed as an upper limit. Therefore we 
observed no clear UV/optical variability accompanying 
the X-ray flare detected in XI. 

3.3. Modeling of the X-ray Spectrum in XI 

We fitted the XI spectrum with several common sim¬ 
ple models. To check the quality of the fits, we re¬ 
binned the spectrum to have a minimum of 20 count 
per bin and adopted the statistic in the fits. How¬ 
ever, all the best-fitting parameters that we will present 
were obtained from the final fits to the spectrum that 
were rebinned to have a minimum of one count per bin 
and adopted the C statistic, considering that the spec¬ 
trum most likely include narrow absorption lines. The 
fits using the C statistic gave smaller uncertainties of 
the parameters than those using the statistic, but 
they are consistent with each other within the uncer¬ 
tainties. Because our source is most likely associated 
with SDSS J152130.72-P074916.5 at z = 0.17901 (Sec¬ 
tion EH), we applied this redshift to the spectral models 
with the convolution model zashift in XSPEC. All mod¬ 
els included the Galactic absorption which w e fixed at 
Vh = 2.61 X 10^° cm“^ (|Kalberla et al.l[200511 using the 
tbabs model. We also included possible absorption intrin¬ 
sic to th e source using the zt babs model. The abundance 
tables of lWilms et al.l (j200r)H were used. 

We first fitted the spectrum with a power law (PL) and 
obtained an unphysically high photon index of PpL = 
5.9±0.2 (the uncertainties from the spectral fits through¬ 
out the paper are all at the 90% confidence level), indi¬ 
cating an ultrasoft source. However, the fit is not accept¬ 
able, with the reduced y^ value of y^ = 1.63 for degrees 
of freedom v = 147. The fits with a MCD model {diskbb 
in XSPEC) is much better, with xt = 1-18. However, 
clear residuals are seen between 0.5-1 keV (middle panel 
in Figured]). We investigated several scenarios to explain 
the residuals. We first tried to add a PL. The fit was only 
marginally improved (at the 90% significance level), if we 
fixed FpL at 2.0 and 3.0, with the PL contributing < 2% 
and < 4% (2cr upper limit) of the rest-frame 0.24-11.8 
keV unabsorbed luminosity, respectively. We next tried 
the emission line explanation for the residuals by adding 
a Gaussian line (gaussian in XSPEC) to the MCD model. 
We obtained a significantly improved fit with xt = 1-09 
for n = 144 (i.e., a total decrease of y^ of 17.8 for three 
more degrees of freedom compared with the MCD fit). 
The final fit using the C statistic is given in Table El Tbe 
line has a centroid rest-frame energy Lga = 0.75lg;Q7 keV 
and a width (jga = 0.02+° °® keV, and the apparent in- 
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ner disk temperature is ^Tmcd = 0.125 ±0.006 keV. The 
neutral absorption is negligible, which is also the case for 
all models tested below. 

We then checked whether the residuals could be due 
to ionized absorption by applying the edge model on the 
MCD model and obtained a good fit with xt = 1-05 
{y = 145), i.e. reduced by 22.5 (or the C statis¬ 
tic reduced by 27.2) for two more degrees of freedom. 
The final fit using the C statistic is given in Table [31 
The edge energy is ifedge = 0.85 ±0.02 keV (rest-frame), 
and optical depth is Tedge = 1.0 ± 0.3. The apparent 
inner disk temperature is now /cTmcd = 0.144 ± 0.007 
keV. To obtain the confidence level of the improvement 
of the fit from adding the edge mo del, we applied the 
posterior predictive p- value method ([Hurkett et al.ll2008l : 
iProtassov et al.l l2002f) . With 10,000 spectra simulated, 
we found none having the reduction of the C statistic 
from the introduction of the edge model to be as large 
as obtained for the XI spectrum (i.e., 27.2). Therefore 
the chance probability of improving the fit by adding the 
edge model is < 10 “^. 

We tried to replace the edge model with the more phys - 
ical warm absorber model zxipcf ([Reeves et al.l 1200811 . 
We obtained a good fit with xt = 1-04 {v = 144). The 
final fit using the C statistic is given in Table [31 The 
warm absorber has IVh = ( 51 * 11 ) ^ 10 ^^ cm“^, the ion¬ 
ization parameter log^ = 2.0 ± 0 . 1 , and an outflow line- 
of-sight velocity of ( 0.12 ± 0 . 02 )c, where c is the speed 
of light. The apparent inner disk temperature is slightly 
higher than inferred from the models investigated above 
(/cTivicD = 0.17 ± 0.02 keV). The fit and the correspond¬ 
ing residuals, using the statistic, are shown in the top 
and bottom panels in Figure SI respectively. For this 
model, we checked the possible presence of a PL com¬ 
ponent and found it to improve the fit only at the 90% 
confidence level and contribute < 4% and < 12% (2cr 
upper limit) of the rest-frame 0.24-11.8 keV unabsorbed 
luminosity, if we fixed Fpl at 2.0 and 3.0, respectively. 

Ultrasoft X-ray spectra are also often fitted with 
a single-temperature blackbody (BB, bbodyrad in 
XSPEC), and our fits with this model are given in Ta¬ 
ble [31 Because the residuals in 0.5-1 keV seen in the 
MCD fit are also present in the BB fit, we also tried a 
Gaussian line, an edge and a warm absorber to account 
for them. Generally, the fits using a BB have a x^ value 
higher by 4-8 than those using a MCD. The BB temper¬ 
atures inferred are in the range between 0.10-0.13 keV. 

Finally, we roughly estimated the BH mass with 
the AGN s pectra l model optxagnf (in XSPEG) by 
[Done et ^ Il2012[ l. assuming that the XI spectrum is 
due to pure thermal disk emission (subject to a warm 
absorber). We note that optxagnf was created assum¬ 
ing a disk inclination of 60°. Table [3] gives the results 
of two fits, one for a non-rotating Schwarzschild BH and 
the other for a maximally rotating Kerr BH. The former 
inferred the BH mass of (1.9 ± 0.4) x 10^ Mq and the 
Eddington ratio of 4.2tio) while the latter inferred the 
mass of (1.4 ± 0.3) x 10® Mq and the Eddington ratio of 
0.51^0 12 - The quality of both fits is as good as that we 
obtained with the MCD model, as indicated by the x^ 
value. 



Fig. 5. — The observer-frame 0.2-2 keV pn background sub¬ 
tracted light curve from XI created with the SAS tool epiclccorr. 
The bin size is 300 s. Data in the strong flare background intervals 
in the middle of the observation are shown in gray and seem to 
smoothly connect the data before and after the background flare. 

3.4. The Short-term Variability in XI 

Figure [5] shows the observer-frame 0.2-2 keV pn light 
curves binned at 300 s from XI. The source is clearly 
brighter at the end of the observation than at the be¬ 
ginning. The probability that the count rate is constant 
(excluding the high background flaring period) is lO”**"^, 
from the x^ test. There may be fast variations by a fac¬ 
tor of ~2 on timescales of thousands of seconds, such as 
those at 2.2 x 10^ s and 4.3 x 10"* s into the observation. 
We extracted a low-state spectrum from the first 28 ks 
of the observation and a high-state spectrum from the 
last 14 ks of the observation, excluding data in the high 
background intervals, and obtained the 0.6-1.2 keV to 
0.2-0.6 keV pn count rate ratio to be 0.157 ± 0.015 and 
0.262 ± 0.021, respectively. Therefore the spectrum is 
softer at lower flux, at the A.la significance level. 

We fitted these low-state and high-state spectra with 
the MCD model subject to a fast-moving warm absorber 
(Section 13.31) . The chance probability of improving the 
fit by adding the warm absorber is still low, 0.0003 and 
0.001 for the low-state and high-state spectra, respec¬ 
tively, based on the posterior predictive p-value method 
using 10,000 simulations for each spectrum. The param¬ 
eters of the neutral and ionized absorption were consis¬ 
tent with and thus fixed at values inferred from the fit 
to the average spectrum of the whole observation. We 
obtained the disk parameters kTyiCD = 0.156 ± 0.008 
keV and Vmcd = 93t^g for the low-state spectrum and 

/cTmcd = 0.189 ± 0.010 keV and Vmcd = dOtJ^ for 
the high-state spectrum. The disk temperatures differ 
by 33 eV at the 4.3 ct confidence level, while the inner 
disk radii differ by 24% at the 2.2a confidence level. The 

mass accretion rates, M oc Tmcd-^mcd MCD 

model, differ by 11%, less than Icr uncertainty (we used 
the XSPEC command steppar on Tmcd and Vmcd to 
calculate the uncertainty of M). Therefore one explana¬ 
tion for the short-term variability in XI is fast change of 
the inner disk radius at an approximately constant mass 
accretion rate into the disk. Figure [5] shows that the 
change occurred within ^^10 ks. 
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4. DISCUSSION AND CONCLUSIONS 

J1521+0749 is an ultrasoft X-ray transient consistent 
with the center of a galaxy at z = 0.17901. This galaxy 
is consistent with an old (> 10 Gyr), passive stellar pop¬ 
ulation, based on our modeling of its SDSS spectrum, 
which shows no clear emission lines. The Scr upper limit 
of the bolometric luminosity of the persistent nuclear ac¬ 
tivity estimated from the [O III] upper limit is nearly 
three orders of magnitude lower than the peak X-ray 
luminosity in XI. Therefore we disfavor AGN activity 
as the explanation for the X-ray flare. Instead, it is a 
good soft TDE candidate. We only have one X-ray de¬ 
tection (XI), which does not allow us to formally test 
whether the event follows a decline. However, we 

have a deep observation G1 relatively close to XI (i.e., 4.7 
months before XI), and the non-detection of the source 
in that observation allows us to constrain the disruption 
time tD to be close to or after Gl. In Figured we plot a 
possible — decay curve, with to assumed to be 

at the time of Cl. The non-detection limits of our source 
in SI and X2, the latter of which is relatively deep, are 
consistent with this curve. 

We fortunately caught the X-ray flare in a deep obser¬ 
vation (XI). We detected short-term flux and spectral 
variability in XI, with the X-ray spectrum being softer at 
lower flux. Based on our spectral modeling, the variabil¬ 
ity could be due to change in the inner disk truncation 
radius, by 24% within ~10 ks, at a constant mass ac¬ 
cretion rate. The variation timescale is thus 40 ks, close 
to the viscosity timescale at the innermost stable circular 
orbit, which is ^70 ks, assuming a viscosity parameter of 
0 . 1 , a ratio of the disk thickness to the inner disk radius 
of 0.1, and a BH mass of 10® Mq (see Equation 5.69 in 
iFrank et al.ll200§l . Therefore the short-term variability 
in XI could be due to viscosity instability of the disk. 
Another soft TDE c andidate 2XMM i J184725.1-631724 
that we discovered in lLin et ahl ()2011h also showed strong 
short-term variability in two deep XMM-Newton obser¬ 
vations in the flare peak. We suggested the cause of 
its short-term variability to be the fast change in the 
mass accretion rate, considering the increase of the in¬ 
ner disk temperature with the flux as inferred from the 
fits to a low-state spectrum and a high-state spectrum 
in the second (brighter) XMM-Newton observation. We 
revisit this problem using the procedure as we adopted 
for J1521-I-0749, by calculating the mass accretion rate 

directly through M ex Tmcd-^mcd MCD when 

fitting the low-state and high-state spectra of the second 
XMM-Newton observation and fixing the absorption col¬ 
umn density at the value inferred from the fit to the av¬ 
erage spectrum from the whole observation Hn lLin et al.1 
(IMl . we allowed the column density to be free but tied 
between the two spectra in fits; fixing the column density 
instead allows to obtain better constraints on the MCD 
parameters and thus more accurate comparison between 
the two spectra). In this way the disk temperatures of 
the low-state and high-state spectra were found to differ 
by 10 eV at the 7.9cr confidence level, the inner disk radii 
by 17% at the 3.3cr confidence level, and the mass accre¬ 
tion rates by 1%, less than la uncertainty. Therefore the 
short-term varialiility in 2XMMi J184725.1-631724 could 
also be due to change in the inner disk truncation radius 
at a constant mass accretion rate caused by disk insta¬ 


bility, as we suggest for J1521-I-0749. The TDE candi¬ 
date SDSS J120136.02-1-300305.5 als o showed short-tem 
variability, but its cause is not clear (|Saxton et al.l[20T^ . 
Galactic stellar-mass BH X-ray binaries in the thermal 
disk domin ated state tend to show very we ak short-term 
variability (|Remillard fc McGlintoclJ 120061 1. The short¬ 
term variability could be a unique phenomenon for soft 
TDEs. 

We also detected some edge residuals in 0.5-1 keV 
when fitting the XI spectrum with a MCD. One expla¬ 
nation for them is the presence of a fast-moving warm 
absorber. This absorber could be the radiation-driven 
outflow expected when the stellar debris fallba cks and ac¬ 
cretes onto the BH at a super-Eddington rate (lReedll98^ 
iStrubbe fc QuataertI l200ii 120111: iLodato fc Rossil 201111 . 
In SDSS .1120136.02-^300305.5. ISaxton et 3 . 11(1201211 also 
detected an edge, at ~0.66 keV, whose cause is not clear. 

We did not detect UV or optical variability between XI 
and X2. Although we expect strong UV/optical emis¬ 
sion from the disk and the outflow in the flare peak, 
such emission could be still too faint to be detected, 
maybe because our OM observations are not very deep 
an d/or the galactic stella r emission is too strong. Based 
on ILodato fc Ros^ (l2011h . the peak (/-band (close to the 
OM B filter) emission has the AB magnitude about 22.3 
mag for a TDE of Mbh within 10®-10® Mq at the dis¬ 
tance of JI52I-I-0749 (i.e., = 866 Mpe), while the 

galaxy has a H-band magnitude of 21.2 ± 0.4 based on 
the X2 observation (Table [2]). The disk UV/optical emis¬ 
sion can also be estimated from modeling the X-ray spec¬ 
trum of XI. For this, it is very important to take into 
account the irradiation of the outer disk by the inner 
disk emission. Therefore we adopted the irradi ated disk 
model diskir in XSPEC (|Gierlihski et al.l|2009ll . We as¬ 
sumed pure thermal disk emission, as the Compton tail 
was found to be very weak. The reddening relation 
E(B — V) = 1.7 X 10“^^ Vh was adopted. The warm 
absorber was also required in the fit, but it was assumed 
to be dust-free. The outer disk radius was set to be 
1000 times of the inner disk radius, which corresponds 
to ~100 times of the tidal disrup tion radius for a BH of 
10® Mq (|Guillochon et all 1201411 . The fraction of bolo¬ 
metric flux thermalized in the outer disk /out is assumed 
to be 0 . 1 , which is unfeasibly large ( /out is typically 
~10~ ® in stellar-mass BH X-ray binaries. iGierlihski et al.l 
1200911 . Even adopting the above assumptions that favor 
the UV/optical emission, we still found that such emis¬ 
sion is about one order of magnitude fainter than the 
galaxy emission (or the detection limit) in the OM filters 
used in XI (Table [5]). 

Because of the super-Eddington accretion in TDEs, rel¬ 
ativistic outflows could be launched and interact with 
interstellar medium, resulting in transient radio syn¬ 
chrotron emission that peaks at '^O.l-l mJy (if the 
sourc e is within Gpe distances) in ^^1 yr after the disrup¬ 
tion ([CiannioF^MetigeBIlQlJ; [Bow3[20Ti|). However, 
despite discovery of two jetted hard TDE candidates, 
there is no clear evidence that soft TDEs also launch 
relativistic jets, because, thus far, there is no radio detec¬ 
tion that can be una mbiguously ascri bed to a relativistic 
jet in soft TDEs fsee lKomossall201^ for a review). We 
searched the Very Large Array (VLA) archival images 
and found only one after 2000 April when our event oc- 
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curred. It was taken on 2000 August 8 (close to XI in 
time) in the L band (1.4 GHz). Our source was not de¬ 
tected, with 3(7 upper limit about 0.8 mJy (three times 
the rms value at the location of our source). The non¬ 
detection is not constraining because the radio emission 
from the jet, even if it was la unched, should be still wea k 
at the beginning of the flare (iGiannios fc Metzgei]l201lll . 

The relatively large positional uncertainty of 
J1521-I-0749 does not rule out that it resulted from 
disruption of a main-sequence star by a massive BH 
residing in an off-nuclear (projected offset <2 kpc) glob¬ 
ular cluster in SDSS J152130.72-1-074916.5. However, 
this explanation should not be very promising because 
if the BH mass is ^10®-10® Mq as inferred from the 
X-ray spectral fit (Section 13.31) the cluster migh t have 
a much higher mass (e.g., iPortegies Zwart et al.l 1200^ 
and thus would be very rare, especially for a host that is 
not very massive. The rate of disrupting a normal star 
by a massive BH in a globular cluster is also predicted 
to be extremely low (<10“^ yr“^ per globular cluster), 
espec ially at late times (>2 Gyr, iBaumgardt et al.l 

HoqI. 

The BH inferred for our source is relatively small but 
might still be a little too large for the white dwarf 
(WD) tidal disruption, which requires the BH mass to 
be <2 X 10® Mq if the disruption needs to be outside 
the event horizon or smaller if the disruption is desired 
to be outside the in nermost stable circular orbit (e.g., 
iKrolik fc Pi^lmTIl . The WD TDE is more prone to 
produce an energetic jet, thus hard X-ray emission, than 
the main-sequence star TDE, because the high density 
of WDs can support a strong magnetic field on the BH 
(|Krolik fc PiranlDOlir i. Even in the case of a slow jet 
when the emission could be dominated by the jet pho- 
tospheric thermal emission, the spectrum is unlikely to 
be as soft as in J1521-I-0749, due to Comptonization 
and b roadening with pair production (iShcherbakov et HI 
120131) . Based on the above considerations, we disfavor 
the WD TDE explanation for J1521-I-0749. 

The supersoft XI spectrum of J1521-I-0749 is hardly 
seen in GRBs, even for the small ultralong class. Only 
the ultralong GRB 060218 showed a similarly so ft spec¬ 
trum (photon index ~5.5. iMargutti et al.l 120151) . The 
main problem of associating J1521-I-0749 with an ultra- 
long GRB is that its host is consistent with an early-type 
galaxy, while most ultralong GRBs have hosts showin g 
intensive star forming activity (e.g., iLevan et al.l 120141) . 
Our source is unlikely to be a SN, because SNe hardly 
have Lx > 10^^ erg s“^ and supersoft X-ray spectra. 
While the shock breakout in SN 2008D had peak Lx 
similar to J1521-I-0749 and show ed relatively soft spectr a 
(maximum photon index ~3.2, iSoderberg et al.l I2008D . 
the event was short (< 600 s) and had a fast rise, ex¬ 
ponential decay profile. In contrast, the XI observation 
of J1521-I-0749 was long (~13 hrs) and had much more 
steady count rates (the second half of the observation had 
0.2-2 keV pn count rates higher by 60% than those in the 
first half only). Although the shock breakout in a massiv e 
star exposion might last for hours (iNakar &: Saril[2M^ . 
the XI observation of J1521-I-0749 was still too long and 
had too steady count rates to be part of such an event; 
the old stellar environment of J1521-I-0749 also argues 
against the explanation of a massive star explosion. 


J1521-I-0749 appears in the direction of the clus¬ 
ter of galaxies MKW 3s {z = 0.045). Its candi¬ 
date host galaxy SDSS J152130.72-1-074916.5 is not a 
member of this cluster. J1521-I-0749 is relatively close 
to, but still too far away to be in, two galaxies in 
the cluster, SDSS J152131.80-1-074851.9 {z = 0.04533 
or cLl = 201 Mpc) and SDSS J152134.15-b074815.3 
(z = 0.04225 or = 187 Mpc, see Figure [p. We 
can rule out J1521-I-0749 as a coronally active star 
based on the ratio of the 0.2-12 keV maximum flux 
to the AT-band flux, which is log(Lx/LiR) = 1.05 for 
J1521-I-0749, much higher than seen in corona lly ac¬ 
tive stars (log(Fx/LjR) < —0.9. iLiii et al.l[2012D . The 
ultrasoft spectrum makes J1521-I-0749 similar to su- 
pers(3ft X-ray source s (SSS; iKahabka k, van den Heuv^ 
[2?M IGreiueH [2?i00 | ). They are mostly (coolin g white 
dwarfs iWDs: iKahabka fc van den Heuvell I2006D or nu¬ 
clear burning of the hydrogen-rich matter on the sur¬ 
face of a WD in the so-c alled close binary supersoft 
sources or supersoft novae (IKahabka fc van den Heuvell 
l2006t iGreined I2000D . J1521-I-0749 is unlikely to be a 
coolin g WD due to its large variability (|Yungelson et al.l 
EMi). It is unlikely to be a close binary supersoft 
source or a supersoft novae either, because such objects 
are uncommon (only a few tens found in our Galaxy, 
IKahabka fc van den Heu'^l2006t lGreined|2000D and the 
chance to find one consistent with the center of a bright 
galaxy is essentially zero. 

Therefore J1521-I-0749 is best explained as a soft TDE. 
Its high-quality X-ray data in XI have allowed us to 
search for short-term variability and carry out detailed 
spectral modeling, which revealed rich information of the 
accretion process in the event. Few TDE candidates have 
such high-quality X-ray data. More soft TDEs of similar 
quality are needed to check whether the short-term vari¬ 
ability and the edge that we detected in J1521-I-0749 are 
common in such events and to understand their origins. 
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